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Abstract. We study the structure of protoneutron stars within the finite-temperature Brueckner-Bethe-Goldstone theoretical 
approach, paying particular attention to how it is joined to a low-density nuclear equation of state. We find a slight sensitivity 
of the minimum value of the protoneutron star mass to the low-density equation of state, whereas the maximum mass is hardly 
affected. 
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1. Introduction 

A lot of effort is currently dedicated to dynamical simulations 
of supernovae explosions (Burrows & Lattimer 1986; Pons et 
al. 1999; Villain et al. 2004; Liebendorfer et al. 2005; Fischer 
et al. 2009) and the subsequent formation and evolution of a 
protoneutron star (PNS). The fundamental input to these calcu- 
lations is the nuclear equation of state (EOS) over a wide range 
of densities, apart from microscopic information regarding dif- 
fusion and cooling processes. The output are time-dependent 
radial profiles of the thermodynamic quantities of interest, such 
as temperature, entropy, particle fractions, etc. 

According to current understanding, one can very roughly 
identify different stages of this process (Burrows & Lattimer 
1986; Pons et al. 1999; Strobel et al. 1999). (i) For about 1 sec- 
ond after supernova core bounce, the system consists of a rela- 
tively cool central region surrounded by a hot mantle, collaps- 
ing and radiating off neutrinos quickly, while still also accret- 
ing material, (ii) For the next 20 or so seconds one can identify 
a slowly developing state of the proper PNS, the system first 
deleptonizing and heating up the interior parts of the star in the 
process, and then beginning to cool down by further neutrino 
diffusion, (iii) After several tens of seconds, the final state of 
a cold NS has basically been formed, which continues to cool 
down slowly first by neutrino and later by photon emission. 

In this article we do not intend to perform dynamical sim- 
ulations, but focus on the consistent construction and evalua- 
tion of basic consequences of the temperature-dependent nu- 
clear EOS during the prominent PNS stage (ii). We therefore 
assume strongly idealized, static profiles representing this en- 
vironment. For simplicity we use a constant entropy per baryon 
throughout the star and investigate the sensitivity of the results 
to the chosen value of entropy S /A, as is usually done (Gondek 
et al. 1997; Goussard et al. 1997). 



Another important aspect is the treatment of neutrino trap- 
ping. Gravitational collapse calculations of the white-dwarf 
core of massive stars indicate that, at the onset of trapping, the 
electron lepton fraction is Y e — x e + x Ve * 0.4, the precise 
value depending on the efficiency of electron capture reactions 
during the initial collapse stage. Also, because no muons are 
present when neutrinos become trapped, the constraint = 
Xfi ~x$ =0 is imposed. We fix the Y\ at these values in our cal- 
culations for neutrino-trapped matter. Furthermore, since low- 
density nuclear matter loses neutrino very fast, the concept of 
a neutrino sphere, inside which neutrinos are trapped, is often 
introduced during the PNS stage (Gondek et al. 1997; Fischer 
et al. 2009). Our treatment of the neutrino sphere is explained 
later in section|4] 

One of the most advanced microscopic approaches to the 
EOS of nuclear matter is the Brueckner theory. In recent years, 
it has made rapid progress in several aspects, (i) The conver- 
gence of the Brueckner-Bethe-Goldstone (BBG) expansion has 
been firmly established (Song et al. 1998; Baldo et al. 2001). 
(ii) The addition of phenomenological three-body forces (TBF) 
based on the Urbana model (Carlson et al. 1983; Schiavilla et 
al. 1986), permitted improving the agreement to a large ex- 
tent with the empirical saturation properties (Baldo et al. 1997; 
Zhou et al. 2004; Li et al. 2006,2008a). (iii) The Brueckner- 
Hartree-Fock (BHF) approach has been extended in a fully 
microscopic and self-consistent way to describe nuclear mat- 
ter also containing hyperons (Schulze et al. 1995, 1998). iv) It 
has also been extended to the finite-temperature regime within 
the Bloch-De Dominicis formalism with a realistic nucleon- 
nucleon interaction (Baldo & Ferreira 1999; Baldo 1999). 

We used the finite-temperature BHF EOS to model PNSs 
in our previous papers (Nicotra et al. 2006; Burgio & Schulze 
2009). The present article represents further evolution in this 
approach, in particular improving two aspects of the previous 
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work: (a) an exact calculation of the finite-temperature EOS, in 
contrast to the so-called "frozen correlations approximation" 
used before; (b) particular attention to how to join the BHF 
high-density EOS to different EOSs describing the low-density 
domain of clustered nuclear matter. 

The structure of this article is as follows. The BHF ap- 
proach at finite temperature is reviewed in Sect. 2, and the 
resulting composition and EOS of beta-stable matter are pre- 
sented in Sect. 3. The joining to different low-density EOSs is 
discussed in Sect. 4, and the final properties of PNSs are shown 
in Sect. 5. Conclusions are drawn in Sect. [6] 

2. Brueckner theory at finite temperature 

The formalism that is closest to the BBG expansion, and that 
actually reduces to it in the zero-temperature limit, is the one 
formulated by Bloch & De Dominicis (1958, 1959a, 1959b). In 
this approach, the essential ingredient is the two-body scatter- 
ing matrix K, which, along with the single-particle potential U, 
satisfies the self-consistent equations 



{k x k 2 \K{W)\hk A ) = {k x k 2 \V\hk A ) 

[l-n(0[l 
+ Re V<*i* 2 |V|*^>- 



n{K)} 



W-Ev- E K + ie 



%k' A \K{W)\hh) 



(1) 



and 



U(k l ) = Y j n(k 2 )(k l k 2 \K(W)\k l k 2 ) A , 



(2) 



where generally denote momentum, spin, and isospin. Here 
V is the two-body interaction, W = + E^ represents the 
starting energy, and Et — k 2 /2m + U(k) the single-particle en- 
ergy. Equation ([TJ coincides with the Brueckner equation for 
the K matrix at zero temperature if the single-particle occupa- 
tion numbers n(k) are taken at T = 0. At finite temperature, 
n(k) is a Fermi distribution. For a given density and tempera- 
ture, Eqs. ([1) and (f2| have to be solved self-consistently along 
with the following equation for the auxiliary chemical potential 



1 



(3) 



The grand-canonical potential density u> in the Bloch-De 
Dominicis framework can be written as the sum of a mean-field 
term and a correlation contribution (Baldo & Ferreira 1999; 
Baldo 1999), 

<a = -V -ln(l + e^-V) + n(k)U(k) 
k \-P 

+ \J ^^~ m K2\nctmi[X(W)n8(Ho-m]j, (4) 

where the trace, Tr2, is taken in the space of antisymmetrized 
two-body states, and the two-body scattering matrix % is de- 
fined by 

{k x k 2 \K(W)\hk A ) = {k x k 2 \K{W)\hk A ) \\ Vl - «/(*) . (5) 

i=l,4 



If one only considers the first term in the expansion of the 
arctan function (Baldo & Ferreira 1999, see Appendix A), then 
the correlation term reduces to 



oj c = Ij]ii(*)I/(*). 



(6) 



which defines the grand-canonical potential in total anal- 
ogy with the BBG binding potential, just using Fermi func- 
tions instead of the usual step functions at zero temperature. 
In this way one neglects a series of terms proportional to 
n(k)[l - n(k)] (or powers of it), which turn out to be negligi- 
ble in the temperature and density ranges relevant for neutron 
and protoneutron stars. 

A further simplification can be achieved by disregarding 
the effects of finite temperature on the single-particle poten- 
tial U, and using the T — results in order to speed up the 
calculations (frozen correlations approximation). This was the 
procedure followed in our previous publications (Nicotra et al. 
2006; Burgio & Schulze 2009). On the contrary, all calcula- 
tions shown in this paper are obtained without introducing any 
approximation on the single-particle potential U, and using the 
exact expression for the grand-canonical potential, Eq. |@}. 

Furthermore, in the previous work the calculations were 
performed using the G-matrix instead of the /if-matrix, [i.e., 
taking the real part in Eq. (f2]i instead of Eq. ([TJ], which pro- 
duces a more repulsive EOS at high density. Strictly speak- 
ing, the Bloch-De Dominicis formalism has been derived for 
the ^-matrix; on the other hand, the use of the G-matrix en- 
sures the compatibility with the usually performed continuous- 
choice BHF calculations at zero temperature. 

In this framework, the free energy density is then 



f = (D+pp.. 



(7) 



and all other thermodynamic quantities of interest can be com- 
puted from it; namely, the "true" chemical potential /i, pressure 
p, entropy density s, and internal energy density e read as 



P 

s 
e 



df 
dp' 
2 d(f/ P ) 
dp 

J± 

dT ' 
f + Ts. 



(8) 

(9) 

(10) 
(11) 



We stress that this procedure permits to fulfill the Hugenholtz- 
Van Hove theorem in the calculation of thermodynamic al quan- 
tities in the Brueckner theory. For an extensive discussion of 
this topic, the reader is referred to Baldo (1999), and references 
therein. 

It is well known that at zero temperature the non-relativistic 
microscopic approaches do not correctly reproduce the nuclear 
matter saturation point po ~ 0.17 frrT 3 , E/A a -16 MeV. One 
common way of correcting this deficiency is to introduce three- 
body forces among nucleons. Since a complete microscopic 
theory of TBF is not available yet, we have adopted the phe- 
nomenological Urbana model (Carlson et al. 1983; Schiavilla 
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Fig. 1. Free energy per nucleon for symmetric (upper panel) 
and pure neutron (lower panel) matter (solid lines). The tem- 
peratures vary from to 50 MeV in steps of 10 MeV. See text 
for details. 



et al. 1986), which consists of an attractive term due to two- 
pion exchange with excitation of an intermediate A resonance 
and a repulsive phenomenological central term. In the BBG ap- 
proach, the TBF is reduced to a density-dependent two-body 
force by averaging over the position of the third particle, as- 
suming that the probability of having two particles at a given 
distance is reduced according to the two-body correlation func- 
tion. The corresponding EOS reproduces the nuclear matter 
saturation point correctly (Baldo et al. 1997; Zhou et al. 2004; 
Li et al. 2008a) and fulfills several requirements from the nu- 
clear phenomenology. In all calculations presented in this paper 
we used the Argonne Vig nucleon-nucleon potential (Wiringa 
et al. 1995) along with the phenomenological Urbana TBF. 

In Fig. [T] we display the free energy vs. the nucleon den- 
sity, obtained following the procedure discussed above, both 
for symmetric and pure neutron matter, and for several values 
of the temperature between and 50 MeV. The dots represent 
the numerical results, and the curves are empirical fits, which 
can be easily implemented in numerical simulations. The thin 
dashed curves represent the calculations obtained with the G- 
matrix and the frozen correlations approximation for the single- 
particle potential U, as adopted in our previous papers (Nicotra 
et al. 2006; Burgio & Schulze 2009). The discrepancy at high 
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Fig. 2. The entropy per particle for symmetric (upper panel) 
and pure neutron (lower panel) matter at different temperatures. 

density is mainly due to the use of the G-matrix instead of the 
Zf-matrix, whereas the differences observed at low density and 
high temperature are consequences of the frozen correlations 
approximation. 

We find that the following functional forms provide excel- 
lent parametrizations of the new numerical results in the re- 
quired ranges of density (0.03 firT 3 < p < 1 firT 3 ) and temper- 
ature (0 MeV < T < 50 MeV) for symmetric (SNM) and pure 
neutron matter (PNM): 

p< 

-(T,p) = -(137 + 157f 2 )p + 308p L82 + 207f 2 ln(p) 

+(-47.5f 2 + 71f 241 )/p-5, (SNM), (12) 

-(T,p) = (11 - 122f 2 )p + 309p 195 + 173f 2 ln(p) 
A 

+(-48f 2 + 71f 2 35 )/p + 6, (PNM), (13) 

where t = T/(100 MeV) and F and p are given in MeV and 
firT 3 , respectively. 

We notice that the free energy of symmetric matter shows 
a typical Van der Waals behavior (with Tq ~ 19 MeV, p c » 
0.06 fm 3 ) and is a monotonically decreasing function of the 
temperature. At T — the free energy coincides with the total 
energy and the corresponding curve is just the usual nuclear 
matter saturation curve. 
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Fig. 3. The free symmetry energy as a function of the nucleon 
density at different temperatures. 



From the free energy we can calculate the entropy (in units 
of the Boltzmann constant) from the thermodynamical relation- 
ship 



(14) 



and this is displayed in Fig.[2]for both SNM and PNM. 

It turns out that the dependence of the free energy on the 
proton fraction can be approximated very well by a quadratic 
dependence, as at zero temperature (Bombaci & Lombardo 
1991; Baldo et al. 1998, 2000a): 

j(r,p,jc)« j(T,p,x = 0.5) + (l -2x) 2 F sym (T,p) , (15) 

where the symmetry energy F sym can be expressed in terms of 
the difference of the energy per particle between pure neutron 
(x = 0) and symmetric (x = 0.5) matter: 



F sym (T,p) 



1 d(F/A) 
~4 dx 



(T,p,0) 



j(T,p,0)-j(T,p,0.5). 



(16) 
(17) 



Therefore, in order to treat the beta-equilibrated case, it is only 
necessary to use the parametrizations Eqs. (1121131 1 of the free 
energy for SNM and PNM. 

For completeness, we display in Fig. [3]the behavior of the 
free symmetry energy for several values of the temperature. 
Full dots represent the numerical values, and the solid curves 
are the results of the polynomial fits. Even at finite tempera- 
ture, the symmetry energy increases nearly linearly with den- 
sity, as is generally the case in the BHF approach (Li et al. 
2006,2008a). However, at low density there are significant de- 
viations due to finite temperature, which might have experi- 
mental implications. The symmetry energy plays a significant 
role not only in the cooling of protoneutron stars, but also in the 
dynamics of heavy ion collisions induced by radioactive beams 
and the structure of exotic nuclei. For a review, see Li et al. 
(2008b). 



3. Composition and EOS of hot stellar matter 

In neutrino-trapped yS-stable nuclear matter, the chemical po- 
tential of any particle i - n,p,l is uniquely determined by the 
conserved quantities baryon number Z?,, electric charge Q/, and 
weak charges (lepton numbers) Lf \ L. : 



pi = Bj/j n - Qi(p„ - p p ) + Lfp Ve + Lf'p 



rW, 



(18) 



For stellar matter containing nucleons and leptons as relevant 
degrees of freedom, the chemical equilibrium conditions read 
explicitly as 



p n - p p = p e - p Ve =Pp_+ Py 



(19) 



At given baryon density p, these equations have to be solved 
with the charge neutrality condition 

2 Q*i = ( 20 ) 

i 

and those expressing conservation of lepton numbers 

Yi = Xi - xj + x Yl - x vt , I — e,p . (21) 

As discussed in Sect. [1] we fix the lepton fractions to Y e = 0.4 
and Yft = for neutrino-trapped matter. When the neutrinos 
have left the system, their partial densities and chemical poten- 
tials vanish and the above equations simplify accordingly. 

The nucleon chemical potentials are obtained from the free 
energy density /, Eq. (0, 



df_ 

dpi 



n, p . 



(22) 



and the chemical potentials of the non-interacting leptons are 
obtained by solving numerically the free Fermi gas model at fi- 
nite temperature. Once the hadronic and leptonic chemical po- 
tentials are known, one can proceed to calculate the composi- 
tion of the /J-stable stellar matter, and then the total pressure p 
through the usual thermodynamical relation 



P=P 



dp 



(23) 



Following this procedure, we first discuss the populations 
of beta-equilibrated stellar matter. In Fig. [4] we display the rel- 
ative particle fractions as a function of the baryon density for 
cold (left panels) and hot beta-stable matter characterized by 
entropy values S /A = 1 (central panels) and S /A = 2 (right 
panels). The upper panels show the particle fractions when 
stellar matter contains only neutrons, protons, electrons, and 
muons, whereas the lower panels show the particle fractions 
in neutrino-trapped matter. We notice that the electron frac- 
tion is greater in neutrino-trapped matter than in the neutrino- 
free case; therefore, the nuclear matter is more symmetric, and 
the resulting EOS will be softer. The appearance of muons in 
trapped matter is shifted to higher density values, because their 
onset is determined by the difference between the neutron and 
proton chemical potentials, Eq. ( fT9] >. 

We observe that thermal effects influence the populations 
mainly in the low-density region. In fact, the presence of tails in 
the Fermi distribution makes it possible to create (anti)particles 
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Fig. 4. Relative populations as functions of baryon density in neutrino-free (upper panels), and neutrino-trapped (lower panels) 
beta-equilibrated matter at entropies S /A = 0, 1, 2. 



at any density and thus typical production thresholds, like for 
muon creation, disappear at finite entropy. At high density, ther- 
mal effects are less important and do not change the composi- 
tion appreciably when increasing entropy. As already found in 
(Nicotra et al. 2006) in the isothermal case and in (Burgio & 
Schulze 2009), the proton fraction increases slightly at high 
density with increasing entropy, thus causing a softening of the 
EOS. This is at variance with the results shown in (Prakash et 
al. 1997), and is presumably caused by the different tempera- 
ture dependence of the potential part of the EOS. 

In Fig. |5j we display the pressure for beta-stable asym- 
metric matter with and without neutrinos, as a function of the 
energy density at entropy S /A = 2 (red curves), in compar- 
ison with cold beta-stable matter (black curves). As discussed 
above, we notice that both thermal effects and neutrino trapping 
produce a softer EOS than in the cold case. 

4. Low-density EOS 

The Brueckner approach provides a realistic modeling of nu- 
clear matter only at densities above about half normal nuclear 



matter density. Below this threshold, clusterization sets in and 
the system becomes inhomogeneous. Another theoretical ap- 
proach therefore has to be used in this "low-density" regime, 
and we employ two widely used liquid-drop-type models at fi- 
nite temperature, namely the one of Lattimer & Swesty (1997) 
and the one of Shen et al. (1998a, 1998b). In the first case 
we in fact compare a model with a (too) large compressibility, 
K = 370 MeV, here denoted by "LS" and a second one with a 
lower compressibility, K = 263 MeV, denoted by "SKa". The 
Shen EOS is characterized by a compressibility K = 28 1 MeV 
and a symmetry energy at saturation of E sym = 36.9 MeV, 
while the BHF EOS has K = 210 MeV and £ svm = 34.7 MeV. 



Of course, since no phase transition is involved, but only 
two different theoretical descriptions of the same state of mat- 
ter, the joining of the two EOSs requires the thermodynami- 
cal observables, i.e., free energy, internal energy, pressure, and 
chemical potentials, to be continuous functions of the baryon 
density. Of particular interest are therefore (i) the consistent 
joining of the low- and high-density EOS, and (ii) eventual 
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Fig. 5. Pressure as a function of energy density for beta- 
equilibrated cold and hot matter with and without neutrino trap- 
ping. 

differences caused by the use of different low-density models. 
These aspects are studied in detail in the following. 

Another important feature of the low-density domain is 
the treatment of neutrino trapping. Physically, neutrinos es- 
cape rapidly from the low-density matter during the PNS evo- 
lution, and so the trapping condition Y e = 0.4 does not apply 
anymore. This effect can be roughly modeled by the concept 
of a neutrino sphere inside which the neutrinos are trapped. 
Typical model-dependent values for the location of the neutrino 
sphere found in the literature are 2 x 10~ 3 for 3 (Gondek et al. 
1997), 6 x 10~ 4 for 3 (Strobel et al. 1999), and 2 x 10 5 for 3 
(Fischer et al. 2009). Given these variations, we choose the fol- 
lowing "natural cutoff" procedure: By imposing the condition 
Y e = 0.4 at any density, at a certain threshold density p v « 
10~ 5 - 10~ 6 fm"" 3 , the electron fraction x e becomes 0.4, and the 
neutrinos disappear naturally. For lower densities we consider 
the matter untrapped. This is illustrated in Fig. [6] showing the 
electron and neutrino fractions following this recipe with the 
different EOSs. This simple procedure avoids making assump- 
tions about the neutrino sphere, but clearly a satisfactory treat- 
ment of neutrino trapping would require coupled dynamical 
calculations of PNS evolution and the temperature-dependent 
EOS, which might be possible in the future. 

In Fig. [7] we display the joint low-density LS, SKa, Shen 
+ BHF EOSs over a wide range of densities. The black lines 
represent the calculations performed at T — 0, whereas the red 
(blue) lines are the ones for fixed entropy S I A = 1 (S /A = 2) in 
neutrino-trapped matter. The transition densities between low- 
density inhomogeneous regime and homogeneous BHF regime 
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Fig. 6. The lepton fractions in beta-stable matter as function of 
baryon density for the different EOSs. 



in the different configurations are displayed in Fig.|7]as vertical 
bars, and summarized in Table Q] One observes in general a re- 
duction of the transition density with decreasing compressibil- 
ity of the low-density EOS and with increasing temperature. In 
all configurations, a consistent joining of the two segments of 
the EOS is possible, fulfilling the above-mentioned continuity 
requirements. In practice we performed a Maxwell construc- 
tion by equating pressure and chemical potentials of the low- 
and high-density sectors, and verified that the other thermody- 
namic variables do not exhibit significant discontinuities at the 
transition point. 

As a further general prominent feature, we notice a strong 
stiffening of the EOS in neutrino-trapped matter as compared to 
the cold untrapped case in the domain p * 10~ 4 , . . . , 10" 1 fm" 3 
(black vs. colored curves in Fig. |7j, which comes from the 
greater lepton Fermi pressure in trapped matter (larger electron 
fraction and additional neutrino contribution, cf., Fig.[6]l. 



Table 1. Transition densities (in fm 
(LS, SKa, Shen) and BHF EOS. 



') between low-density 





untrapped, T=0 


trapped, S/A=l 


trapped, S/A=2 


LS 


0.160 


0.145 


0.130 


SKa 


0.110 


0.090 


0.086 


Shen 


0.120 


0.110 


0.094 
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Fig. 7. The EOS for cold untrapped matter (black curves) and 
for neutrino-trapped matter with entropy S /A = 1 (red curves) 
and S /A =2 (blue curves). 

5. (Proto)neutron star configurations 

The stable configurations of a (P)NS can be obtained from 
the well-known hydrostatic equilibrium equations of Tolman, 
Oppenheimer, and Volkov (Shapiro & Teukolsky 1983) for 
pressure p(r) and enclosed mass m(r), 

dp __ Gme ( ! + Pl e ){ 1 + ^p/m) 
dr r 2 1 - 2Gm/r 

dm 
dr 



Anr 2 e . 



(24) 



(25) 



once the EOS p{e) is specified, with e — e^+ei the total internal 
energy density (G is the gravitational constant). For a given 
central value of the energy density, the numerical integration 
of Eqs. ( l24b and d25l l provides the mass-radius relation. 

Our results, using the different EOSs introduced in the pre- 
vious section, are displayed in Fig. [8] which shows the gravi- 
tational mass (in units of the solar mass M = 1.98 x 10 33 g) 
as a function of the radius (left panels), and the central density 
(right panels) (normalized with respect to the saturation value 
eo = 152 MeV fm~ 3 ). The black curves represent the calcula- 
tions performed at T — 0, whereas the thick (thin) colored lines 
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Fig. 8. The gravitational mass as a function of the radius (left 
panels), and the central energy density (right panels), obtained 
with the LS (top panels), SKa (central panels), and Shen (bot- 
tom panels) EOSs. 



denote stable configurations in neutrino-free (neutrino-trapped) 
matter at constant entropy S /A — 1 (dashed lines) and S /A — 2 
(dot-dashed lines). 

We notice that the value of the maximum mass decreases 
slightly in neutrino-trapped matter, because of the softening of 
the EOS discussed before. The value of the maximum mass 
turns out to be practically independent of the entropy, whereas 
the value of the radius for a fixed gravitational mass depends 
strongly on the entropy and on the presence of neutrinos. The 
dots mark the values of the minimum mass and the correspond- 
ing radius a PNS may have, with typical values in the range 
0.6-0.77 M and about 40-50 km, according to the chosen low- 
density EOS. The proper values are summarized in Table [2] 
One notes that the minimum mass configurations only probe 
the low-density part of the EOS. As expected, the minimum 
mass increases with decreasing compressibility of the EOS, 
and with increasing entropy of the matter. 
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Table 2. Properties of (P)NS minimum and maximum mass 
configurations. 







minimum mass 


maximum mass 






M/M Q 


R(km) 


Pel Pl) 


M/M e 


R(km) 


Pc/PO 


untrapped 


LS 








2.03 


9.86 


10.55 


T=0 


SKa 








2.03 


9.86 


10.42 




Shen 








2.03 


9.93 


10.42 


trapped 


LS 


0.58 


40 


1.02 


1.95 


10.2 


11.34 


S/A=l 


SKa 


0.60 


38 


1.08 


1.95 


10.2 


11.20 




Shen 


0.58 


44 


1.02 


1.95 


10.3 


11.20 


trapped 


LS 


0.70 


44 


0.90 


1.95 


10.7 


10.85 


S/A=2 


SKa 


0.77 


42 


0.90 


1.95 


10.8 


10.70 




Shen 


0.75 


52 


0.77 


1.95 


10.8 


10.80 



Our values of the minimum mass are very similar to the 
ones found in Gondek et al. (1997), where the LS EOS was 
used over the full density range of the PNS. On the other hand, 
our values turn out to be lower than those calculated in Strobel 
et al. (1999), the reason being that our minimum mass is calcu- 
lated during stage (ii) in the PNS evolution, when neutrinos are 
still trapped, and the core and the mantle are characterized by 
the same entropy values as described in the Introduction. 

6. Conclusions 

We have presented a microscopic calculation of the EOS of hot 
nuclear matter in the BHF approach and provided convenient 
parametrizations of the free energy as a function of baryon den- 
sity, proton fraction, and temperature. All other thermodynamic 
quantities of interest can be derived from these. This EOS was 
joined to different standard EOSs describing the low-density, 
inhomogeneous domain of hot nuclear matter, and we roughly 
modeled the structure of PNSs, using strongly idealized tem- 
perature and neutrino-trapping profiles. 

Altogether we find only small variations in our results for 
the minimum mass configurations when different low-density 
EOSs are employed, while maximum masses are practically 
unaffected. Both neutrino trapping and (to a minor degree) fi- 
nite temperature decrease the asymmetry of nuclear matter, and 
thus soften the EOS and decrease the maximum mass. 

However, to provide more quantitative results, it will be 
necessary to perform consistent dynamical simulations of the 
PNS evolution, using self-consistent temperature and neutrino- 
trapping profiles, along with the temperature dependent EOS. 
This is a difficult task for the future. 
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